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Vortex tubes in turbulence velocity fields at Reynolds numbers ReA — 300—1300 
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The most elementary structures of turbulence, i.e., vortex tubes, are studied using velocity data 
obtained in a laboratory experiment for boundary layers with Reynolds numbers Re^ = 295-1258. 
We conduct conditional averaging for enhancements of a small-scale velocity increment and obtain 
the typical velocity profile for vortex tubes. Their radii are of the order of the Kolmogorov length. 
Their circulation velocities are of the order of the root-mean-square velocity fluctuation. We also 
obtain the distribution of the interval between successive enhancements of the velocity increment as 
the measure of the spatial distribution of vortex tubes. They tend to cluster together below about 
the integral length and more significantly below about the Taylor microscale. These properties are 
independent of the Reynolds number and are hence expected to be universal. 

PACS numbers: 47.27.Ak, 47.27.Jv, 47.27.Nz 



I. INTRODUCTION 

Turbulence contains vortex tubes as the most elemen- 
tary structures 0, S S H- Regions of strong vortic- 
ity are organized into tubes. The energy dissipation 
is significant around them. They occupy a small frac- 
tion of the volume and are embedded in the background 
flow that is random and of large scales. Their existence 
has been established at microscale Reynolds numbers 
Re A ^ 2000, by seeding a turbulent liquid with gas bub- 
bles and thereby visualizing regions of low pressure that 
are associated with strong vorticity [E IE 13 • 

Direct numerical simulations 0, |E Hol Hlf have de- 
rived basic parameters of vortex tubes at low Reynolds 
numbers, Rca ^ 200. The radii are of the order of the 
Kolmogorov length rj. The total lengths are of the or- 
der of the integral length L. The circulation velocities 
are of the order of the root-mean-square velocity fluctu- 
ation {u^y^"^ or the Kolmogorov velocity uk- Here (•) 
denotes an average. The lifetimes are of the order of the 
large-eddy turnover time L/(u^)^/^. 

The universality of these tube parameters has not been 
established because their behavior at high Reynolds num- 
bers has not been known. It is difflcult to conduct a direct 
numerical simulation at ReA ^ 200 il^]. Bubble visual- 
ization does not provide sufficient information jl3j . We 
accordingly use velocity fields of laboratory turbulence at 
high Reynolds numbers to study some of the tube param- 
eters. The velocity field is intermittent at small scales. A 
small-scale velocity variation is enhanced at the positions 
of vortex tubes |E[llL5.I6.,i7, 18, 19]. 

There are several possible configurations for laboratory 
experiments. Although the most popular configuration 
is isotropic grid turbulence, its Reynolds number is not 
high 131 . We instead use rough- wall boundary- layer tur- 
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bulence. The highest Reynolds number achieved in our 
experiment is Rca = 1258, which exceeds those in almost 
all t he p revious studies for vortex tubes in the velocity 
field m. 

Using a probe suspended in the flow, we obtained a 
one-dimensional cut of the velocity field. We measured 
not only the velocity component in the mean-flow direc- 
tion but also the component that is perpendicular to the 
mean-flow direction. The latter component is suited to 
detecting circulation flows such as those associated with 
vortex tubes CI El El E [13 . 

The experiment is described in Sec. We present 
a model for vortex tubes in Sec. IIIII From the experi- 
mental data, the typical velocity profile for vortex tubes 
is extracted and its radius and circulation velocity are 
obtained in Sec. IIVI The spatial distribution of vortex 
tubes is obtained in Sec. The dependences of these 
tube parameters on the Reynolds number are studied in 
Sec. IVil We conclude with remarks in Sec. IV 111 



II. EXPERIMENT 

The experiment was done in a wind tunnel of the Mete- 
orological Research Institute. We use the coordinates x, 
y, and z in the streamwise, spanwise, and floor-normal di- 
rections. The corresponding wind velocities are u, v, and 
w. The origin x = y = z = Ois taken on the tunnel floor 
at the entrance to the test section. Its size was 5x = 18 
m, (5y = 3 m, and Sz = 2 m. Over the entire floor of 
the test section, we placed blocks as roughness elements. 
Their size was 6x = 6 cm, Sy = 21 cm, and Sz = 11 cm. 
The spacing of adjacent blocks was Sx = 6y = 0.5 m. We 
set the incoming-wind velocity to be Ui = 2, 4, 8, 12, 16, 
or 20 m s~^. 

The streamwise and spanwise velocities were simulta- 
neously measured using a hot-wire anemometer. The 
anemometer was composed of a crossed-wire probe and 
a constant temperature system. The wires were made of 
platinum-coated tungsten, 5 /xm in diameter, 1.25 mm in 
effective length, 1 mm in separation, 280 °C in tempera- 
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TABLE I; Summary of experimental conditions, flow characteristics, and parameters of vortex tubes. The kinematic viscosity 
reflects the air temperature in the wind tunnel. The velocity derivative was obtained as — [8v{x + 5x) — 8v{x — Sx) — v{x + 
2Sx) + v{x - 2Sx)]/l2Sx with Sx = U/fs- 



quantity unit values for individual data sets 



experimental conditions: 


















incoming-wind velocity 


Ui 


m s ^ 


2 


4 


8 


12 


16 


20 


sampling frequency 


fs 


kHz 


4 


8 


16 


24 


32 


40 


kinematic viscosity 


V 




0.144 


0.145 


0.145 


0.147 


0.150 


0.149 


99% thickness 




cm 


78 


80 


80 


79 


79 


79 


displacement thickness 


J^f{l - U/tj)dz 


cm 


20 


23 


23 


23 


22 


22 


measurement height 




cm 


35 


35 


30 


30 


25 


25 


flow characteristics: 


















mean streamwise velocity 


U 


m ^ 


1..59 


3.08 


5.83 


8.81 


11.1 


13.8 


streamwise flatness factor 






2.73 


2.68 


2.70 


2.69 


2.70 


2.71 


spanwise flatness factor 






3.03 


3.03 


3.00 


3.06 


3.01 


3.02 


energy dissipation rate 


(e) = 15i^{{a^vf)/2 


2 -3 

m s 


0.0316 


0.254 


2.03 


5.65 


14.6 


26.2 


streamwise velocity fluctuation 




m s~ ^ 


0.271 


0.554 


1.15 


1.71 


2.38 


2.95 


spanwise velocity fluctuation 




m s~ ^ 


0.227 


0.462 


0.958 


1.42 


2.01 


2.53 


Kolmogorov velocity 


UK = (i^ie))^^^ 


m s~ ^ 


0.0260 


0.0438 


0.0737 


0.0955 


0.122 


0.141 


streamwise integral length 


Lu — f {u{x + 5x)u{x)) I {u^)d5x 


cm 


48.4 


44.2 


39.0 


42.3 


43.6 


43.6 


spanwise integral length 


= J {v(x + Sx)v(x))/{v^)dSx 


cm 


7.76 


7.08 


7.06 


6.23 


5.79 


5.92 


Taylor microscalc 




cm 


1.88 


1.35 


0.991 


0.889 


0.788 


0.740 


Kolmogorov length 




cm 


0.0554 


0.0331 


0.0197 


0.0154 


0.0123 


0.0106 


microscalc Reynolds number 






295 


430 


655 


861 


1054 


1258 


parameters of vortex tubes: 


















radius 


Ro 




6.08 


6.04 


6.28 


7.14 


6.98 


7.37 


circulation velocity 


Vo 




0.600 


0.526 


0.485 


0.476 


0.469 


0.464 


circulation velocity 


Vo 


UK 


5.23 


5.55 


6.31 


7.10 


7.73 


8.37 


Reynolds number 


Ron = VoRo/f 




32 


34 


40 


51 


54 


62 


clustering scale 


Sxa 


X 


2.38 


2.90 


2.82 


2.77 


2.67 


2.92 


probability density 






0.349 


0.451 


0.536 


0.508 


0.523 


0.545 



ture, and oriented at ±45° to the streamwise direction. 
The calibration was done before and after each of the 
measurements. We did not measure the floor-normal ve- 
locity, which suffers from mean shear and hence is less 
suited to studying vortex tubes than the spanwise veloc- 
ity. 

The measurement positions were at a; = 12.5 m, where 
the boundary layer was well developed. The 99% thick- 
ness, i.e., the height at which the mean streamwise ve- 
locity U is 99% of its maximum value U, was 0.8 m. The 
displacement thickness /q^(1 — U/U)dz was 0.2 m. Here 

z is the height for the velocity U The 99% thick- 

ness and displacement thickness were independent of the 
incoming-wind velocity. Thus, among all the measure- 
ments, the overall flow structure was the same. The rea- 
son is that the Reynolds number for the entire boundary 
layer was high enough. 

We determined the measurement height z so that the 
flatness factor for the spanwise velocity {v'^) / {v^)"^ was 
close to the Gaussian value of 3. The flatness factor 
was less than 3 at small heights where the flow was af- 
fected by the surface roughness. The flatness factor was 
greater than 3 at large heights where the flow was af- 
fected by the fluctuation of the interface to the outer 
laminar flow. We obtained the Gaussian value of 3 at an 



intermediate height, where eddies with various sizes and 
strengths filled the space randomly and independently 
[2l| . There the flatness factor for the streamwise veloc- 
ity (u^) /(it^)^ was different from 3 because the turbulence 
was not isotropic at large scales. 

The signal was low-pass filtered at /c = 2-20 kHz with 
24 dB per octave and sampled digitally at fs = 4-40 kHz 
with 16-bit resolution. To avoid aliasing, the sampling 
frequency was set to be twice the filter cutoff frequency, 
Js — "^fc- The data length was 2 x 10^ points for the 
incoming- wind velocities Ui = 4, 8, and 12 m s~^. It was 
8 X 10^ points for the incoming-wind velocities Ui — 2, 
16, and 20 m s^^. 

The energy spectra of the spanwise velocity are shown 
in Fig. n We have used Taylor's frozen-eddy hypothesis 
to convert temporal variations into spatial variations in 
the streamwise direction. Throughout the energy spec- 
tra, the signal-to-noise ratio is high. 

The experimental conditions and flow characteristics 
are summarized in Table ^ The microscalc Reynolds 
number Re^ ranges from 295 to 1258. We have ob- 
tained the smallest-scale statistics from the spanwise- 
velocity gradient d^v instead of the usual streamwise- 
velocity gradient d^u, by assuming the smallest-scale 
isotropy {{dxvY) = 2{{dxu)'^). This is because, espe- 
cially at small scales in strong turbulence, the u compo- 
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FIG. 1: Energy spectrum of the spanwise velocity at Re^ = 
295, 430, 655, 861, 1054, and 1258 (from bottom to top). The 
wave number k is in units of m~^ instead of the usual radian 
m~^. The dotted line denotes Kolmogorov's k~^^^ law for the 
inertial range. 



nent measured by a crossed-wire probe is contaminated 
with the w component that is perpendicular to the plane 
of the two wires of the probe 22, 23] . The v component 
is free from such contamination. For the smallest-scale 
isotropy, we have no direct evidence. The smallest-scale 
isotropy still serves as a meaningful assumption, even if 
it was not achieved, because our present results are then 
comparable with those obtained in isotropic turbulence. 



of the circulation flow uq along the x axis are 

A cos 6* 



u{x) 



R 



-ue{R), 



xcosU 
v{x} = ^ U0{R}, 



(2a) 
(2b) 



with 



i?2 ^ x2(l-sin2 6lcosV)-fA2(l-sin2 6lsinV) 

-|-2a;A sin^ 6* sin (yscos (3) 

For the radial inflow uji of the strain field, the streamwise 
and spanwise components are 



u{x) 



v{x) 



x{l — sin^ 9 cos^ ip) + A sin^ 9 sin ip cos (fi 
R 

X siri^ 9 sin ip cos ip + A{1 — sin^ 9 sin^ ip) 
R 



ur{R), (4a) 



uniR). (4b) 



If the vortex tube passes close to the probe (A < Rq) 
and the tube is not heavily inclined (0 ~ 0), the spanwise 
velocity is dominated by the small-scale circulation flow 
[Eq. Ij2bp ] . The streamwise velocity is dominated by the 
large-scale radial inflow [Eq. (|4a|l ] . This situation is of 
our interest. The velocity profiles of vortex tubes with 
A < i?o and ~ are nearly the same Ti]. If A > i?o 
OT 9 ^ 0, the tube signal is weak at least in the spanwise 
velocity at small scales. 

Since there are large-scale velocity fluctuations, a vor- 
tex tube in actual turbulence does not necessarily pass 
the probe along the mean streamwise direction . This 
fact was nevertheless not serious in our experiment. The 
incident angle 4> of the vortex tube was not large, i.e., 
((/)) ~ arctan((w2)i/2/[/) < 10° for (vY^^/U < 0.2 as in 
Table n 



IV. VELOCITY PROFILE OF VORTEX TUBES 



III. MODEL FOR VORTEX TUBES 

The representative model for vortex tubes is the Burg- 
ers vortex, an axisymmetric steady circulation in a strain 
field. In cylindrical coordinates, the circulation and the 
strain field are written, respectively, as 



UQ cx 



aR 



1 — exp 



oR^ 



{ur, ue, uz) = ( ^2"^' 0> ) • 



(a>0), (la) 
(lb) 



Here v is the kinematic viscosity. The circulation is max- 
imal at i? = i?o = 2.24(z^/a)^/2. Thus i?o is regarded as 
the tube radius. 

Suppose that the axis of the vortex tube penetrates the 
{x,y) plane at the point (0, A). The x and y axes are, 
respectively, in the streamwise and spanwise directions. 
If the direction of the tube axis is (0, 1^9) in spherical coor- 
dinates, the streamwise (w) and spanwise {v) components 



The typical profiles for vortex tubes in the streamwise 
(m) and spanwise (w) velocities are extracted by averaging 
signals centered at the position where the absolute value 
of the spanwise- velocity increment \v{x-\-6x) — v(x) \ is en- 
hanced above a certain threshold 0, UtI ITq . We set the 
scale 5x to be the sampling interval U / fs- The threshold 
is set to be the highest percentile for the absolute values 
of the velocity increments. Thus 1% of them are used for 
the averaging. When the velocity increment is negative, 
we invert the sign of the v signal before the averaging. 
The results are shown in Fig. [3 

The threshold value for the enhancement of the veloc- 
ity increment has been determined with a compromise. 
If the threshold is higher, the statistical uncertainty is 
more significant. If the threshold is lower, the contami- 
nation with the background flow is more significant. We 
nevertheless expect that our following results are qual- 
itatively independent of the threshold if the fraction of 
the velocity increments used for the averaging is ^ 1%. 
They comprise the tail of the probability density distri- 
bution that is well above the Gaussian distribution with 
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FIG. 2: Typical profiles for vortex tubes in the streamwise 
(m) and spanwise (v) velocities at Re^ — 295, 430, 655, 861, 
1054, and 1258 (from top to bottom). The u profile is shown 
separately for u(2; + C///s) —M(a;) > and u(a:; + L''//s) — ti(a::) < 
at a; = (designated as u"*" and u~). The position x is 
normalized by the Kolmogorov length rj. On the ordinate, the 
unit scale corresponds to one-tenth of the root-mean-square 
velocity fluctuation (v'^)^^^ . The v profiles of Burgers vortices 
with A = and ^ = are shown with dotted lines [see Eqs. 
ll2bl and ©]. 



FIG. 3: Probability density distribution of the absolute ve- 
locity increment \v{x + U/fa) — v(x)\ at Re^ = 295, 430, 
655, 861, 1054, and 1258 (from top to bottom). The dis- 
tribution is shifted vertically by a factor 10"^. The velocity 
increment is normalized by its standard deviation, which is 
0.00678, 0.0186, 0.0497, 0.0830, 0.124, and 0.167 m s'^ (from 
top to bottom). The arrows indicate the range of the en- 
hanced velocity increments used in our analyses, which share 
1% of the total. The dotted lines denote a Gaussian distribu- 
tion. 



the same standard deviation as shown in Fig. O "24] . The 
only deficit is that the threshold is too high for some weak 
vortex tubes. They are not considered here. 

The u profiles in Fig. |21are separated for u{x + Sx) — 
u{x) > ein.du{x + dx) — u{x) < at a; = (designated as 
w+ and u~). We have decomposed the individual profiles 
into the symmetric and antisymmetric components and 
have shown only the antisymmetric components. The 
contamination with the w component leads to a sym- 
metric positive excursion in the u profile '23', 25|. Such 
an excursion was seen in our previous study ,18, ..26] . The 
u profile averaged for vortex tubes is antisymmetric (see 
Sec.lnl. 

To ensure statistical significance, we have reduced the 



sample size by a factor of 2 and redone the conditional 
averaging. The velocity profiles do not differ from those 
in Fig. lUby more than a line thickness. 

For reference, we show the v profile (j2b|) of a Burgers 
vortex with A = and 6* = in Fig. |21(dotted lines). The 
radius Rq and the maximum circulation velocity Vq have 
been determined so as to fit the observed profile around 
its peaks. There our conditional averaging prefers vortex 
tubes with A ~ and 61 ~ because such tubes have 
the strongest signals. Since the Kolmogorov length rj is 
smaller than the probe size in the streamwise direction 
I, we have assumed that the measured velocity Vm.(x) is 
the true velocity vt (x) averaged over the probe size |22| : 
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Vmix) 



1/2 



Vt{x + x')dx' with I 



-1/2 



1.25 mm 



(5) 



The i?o and T^) values are shown in Tabled While the 
radius i?o is several times the Kolmogorov length ?], the 
circulation velocity Vq is about a half of the root-mean- 



square velocity fluctuation 
Kolmogorov velocity uk 0, 



"^yj^ or several times the 

Eamiiiiiieiiiiii. 

The observed v profile is close to the profile of a Burg- 
ers vortex [nl, [13, [13 ■ We thus confirm the existence of 
vortex tubes and their responsibility for small-scale in- 
termittency at high microscale Reynolds numbers. The 
observed v profile has more pronounced tails than the 
profile of a Burgers vortex. There should be a contribu- 
tion from vortex tubes that are heavily inclined to the 
streamwise direction with p^. 

The profiles are dominated by the circulation fiow 
We of vortex tubes passing the probe at some distances 
A > [Eq. (|2al) ] or with some incident angles </> > 
[T^ . If the radial inflow uji of the strain field were pre- 
dominant, we would observe u(x) ~ —ax/2 [Eq. H4a|l ]. 
This is not the case. The u~ profile only has a somewhat 
larger amplitude than the profile. Unlike a Burgers 
vortex, an actual vortex tube is not always oriented in 
the stretching direction 0, IHIll El 113 • 



SPATIAL DISTRIBUTION OF VORTEX 
TUBES 



J2 
CO 

O 




10 15 20 25 
interval 



Taylor microscale 



The spatial distribution of vortex tubes on a one- 
dimensional cut of a turbulent flow is studied using the 
probability density Pq of the interval 5x' between suc- 
cessive enhancements of the spanwise- velocity increment. 
The enhancement is defined with the same threshold as 
in Sec. II VI We show the probability density distribution 
in Fig. ^ where the interval is normalized by the Taylor 
microscale A in order to cover both the small and large 
intervals. The statistical significance is less than that of 
the velocity pro files in Fig. [21 but it is still sufficient for 
our analysis [23. 

For intervals 5x' ^ 5A, we successfully model the prob- 
ability density with a sum of two exponential functions 
(dotted curves): 

/ 5x' \ f 5x'\ 

Folf^a;') = coexp ( - — j +ciexp ( - — j , (6) 

with 5xq < 6x1. The second term implies that the prob- 
ability density distribution has an exponential tail that 
appears linear on a semi-logarithmic plot (dotted straight 
lines). This is characteristic of the Poisson process of ran- 
dom and independent events [2^. Thus the large-scale 
distribution of vortex tubes is random and independent. 
The first term implies that, with decreasing interval, the 
probability density becomes enhanced over that for the 



FIG. 4: Probability density distribution of the interval Sx' 
between vortex tubes at ReA = 295, 430, 655, 861, 1054, and 
1258 (from top to bottom). The probability density distri- 
bution is normalized by the amplitude of its exponential tail 
(dotted straight lines), and it is shifted vertically by a fac- 
tor 10. The interval is normalized by the Taylor microscale 
A. The dotted curves denote the results of a least-squares fit 
with a sum of two exponential functions at 5x' — 5A-25A [Eq. 
@]. The arrows denote the spanwise integral length Lv 



exponential distribution H H IT^ IT^ . Since the en- 
hancement is significant below about the spanwise inte- 
gral length L„, this is attributable to clustering of vortex 
tubes below the energy-containing scale. It was actually 
demonstrated in direct numerical simulations that strong 
vortex tubes lie on the borders of energy-containing ed- 
dies The clustering scale Sxq is a few times the 
Taylor microscale A (Table QJ. 

With decreasing interval 6x' below about the scale Sxq, 
the probability density becomes enhanced over our model 
(0. The clustering of vortex tubes becomes significant. 
In Table m we show the probability density at 6x' = A, 
which has been normalized by the amplitude of the ex- 
ponential tail. The probability density at the smaller in- 
tervals is not so useful because a very strong vortex tube 
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FIG. 5: Dependence of tube parameters on Re^. (a) Ro/rj. 
{h)Vo/{vy^''. {c)Vo/uK. (d)Reo/Rey^ (e) Reo. (f) fo(A). 
They are individually normalized by the values at Re a = 430. 
The filled circles denote the data at Re^ = 295, 430, 655, 
861, 1054, and 1258 from our present experiment. The open 
circles denote the data at Rca — 105, 165, 225, 292, and 329 
from our previous experiment of grid turbulence p^ . 

could cause more than one enhancement of the velocity 
increment. 



VI. DEPENDENCE ON REYNOLDS NUMBER 

Thus far we have obtained the parameters of vortex 
tubes, i.e., the radius Rq, the maximum circulation ve- 
locity Vb, and the interval distribution Pq. Their de- 
pendences on the microscale Reynolds number Hex are 
studied here. To extend the Ke\ range, we also use veloc- 
ity data from our previous experiment of grid turbulence 
at Hex = 105-329 ^3 • These data are reanalyzed in the 
same manner as for our present data. The results are 
summarized in Fig. [S] where quantities are normalized 
by their values at Re\ ~ 430. 

The tube radius Rq scales with the Kolmogorov length 
r] as Rq (X r] over the entire range of the Reynolds number 
[Fig. EJa)]. This is the most significant scaling law among 
those studied here. 

The circulation velocity Vq scales with the root-mean- 
square velocity fluctuation (ii^)^/^ as Vb cx (w^)^/^ at 
R-Ga ^ 400 [Fig. Elb)]. Since this is not the case at 



Re A ^ 400, the scaling is achieved at high Reynolds 
numbers. Although the velocity fluctuation is a quantity 
for large scales, vortex tubes could be formed via shear 
instabilities on the borders of energy-containing eddies 
0, 0, 01 , where a velocity variation over a small scale 
such as the tube radius could be comparable to the ve- 
locity fluctuation. The circulation velocity also scales 
with the Kolmogorov velocity uk as Vb oc uk [Fig. EIc)]. 
However, at Rca ^ 400, this scaling is less significant 
than the scaling with the velocity fluctuation. 

The scaling laws for the radius Rq and the circulation 
velocity Vb lead to the scaling law for the Reynolds num- 
ber Rcq — RqVq/v that characterizes the circulations of 
vortex tubes: 

Reo (X ReJ^''^ if i?o oc 77 and Vb oc {v^)^^^ , (7a) 
Reo = constant if Ro oc rj and Vo oc uk- (7b) 

These relations are from the definitions of the Reynolds 
number Rca, the Kolmogorov length 77, and the Kol- 
mogorov velocity uk as in Table ^ The present result 
favors the former scaling [Fig. [S{d)] rather than the lat- 
ter [Fig. Ele)] at least for Rca ^ 400. With an increase 
of the Reynolds number Rca, vortex tubes have higher 
Reynolds numbers Rcq and are accordingly more unsta- 
ble |10|| . They would nevertheless survive long enough 
to be observable as distinct entities that are responsible 
for small-scale intermittency. Turbulence is known to be 
more intermittent at a higher Reynolds number Rca 0- 
For general vortex tubes, we do not necessarily expect 

the scaling laws Vq oc {v^y^^ and Rcq oc Rca^^. Weak 
vortex tubes are not considered here because our veloc- 
ity profiles were obtained for enhancements of a velocity 
increment. Actually in direct numerical simulations, the 
scaling law Vq oc uk was obtained when vortex tubes 
were identified as local minima of the pressure 3. The 
scaling law Vb oc (w^)^/^ was obtained when vortex tubes 
were identified as enhancements of the vorticity above a 
threshold T^. 

The probability density Pq at the interval dx' = A ap- 
pears to increase as the Reynolds number increases up 
to Rca ^ 400 [Fig. EJf)]. Above the Reynolds number 
Rca — 400, the probability density is constant. The clus- 
tering of vortex tubes appears to become significant and 
then attain an asymptotic state. There is a discontinuity 
between grid turbulence at Rca ^ 300 and boundary- 
layer turbulence at Re a ^ 300. The spatial distribu- 
tion of vortex tubes is affected by large-scale structures 
of turbulence, i.e., energy-containing eddies. For large 
intervals Sx' ^ A, the probability density distributions 
collapse to single curves according to the experimental 
configurations as shown in Fig. El 

VII. CONCLUDING REMARKS 

The streamwise (u) and spanwise (v) velocities were 
measured simultaneously in rough-wall boundary layers 
with microscale Reynolds numbers Rca — 295-1258. We 
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FIG. 6: Probability density distribution of the interval Sx' 
between vortex tubes in grid turbulence at ReA = 105 ob- 
tained in a 3 X 0.8 x 0.8 wind tunnel, grid turbulence at 
ReA = 165, 225, 292, and 329 obtained in a 18 x 3 x 2 
wind tunnel, and boundary-layer turbulence at ReA ~ 295, 
430, 655, 861, 1054, and 1258. We normalize the probabil- 
ity density distribution by the amplitude of its exponential 
tail. The interval is normalized by the Taylor microscale A. 
The data of boundary-layer turbulence are from our present 
experiment. The data of grid turbulence are from our pre- 
vious experiment Among the data of boundary-layer 
turbulence, those at ReA = 295 yield the lowest probability 
density. 

have used the velocity data to study vortex tubes, the 
most elementary structures of turbulence. 

We have extracted the typical v profile for vortex tubes 



(Fig. 12) . The profile is close to the velocity profile of a 
Burgers vortex. The radius Rq is several times the Kol- 
mogorov length 77. The maximum circulation velocity Vq 
is about a half of the root-mean-square velocity fluctua- 
tion {v^y^^. We have also studied the probability density 
distribution of the interval between vortex tubes (Fig.^. 
The probability density is enhanced below about the inte- 
gral length and more significantly below about the Taylor 
microscale, reflecting clustering of vortex tubes. 

For the first time at high Reynolds numbers, we have 
obtained the scahng laws Rq oc r], Vq (x (v^)^/^, and 

1 /2 

Rcq ~ VoRo/iy cx Re^' . The small-scale spatial distri- 
bution of vortex tubes is the same (Fig. |3J|. Since these 
properties do not necessarily exist at Re^ ^ 400, they 
are achieved asymptotically at Re^ ^ 400. They are ex- 
pected to be universal among vortex tubes in turbulence 
at high Reynolds numbers. To confirm this expectation, 
experiments at the higher Reynolds numbers are desir- 
able. Those at similar Reynolds numbers but under dif- 
ferent experimental configurations are also desirable. 

The vortex tubes have been identified using enhance- 
ments of a velocity increment above a threshold. Thus 
our results are biased against weak tubes. The develop- 
ment of a method to identify vortex tubes with various 
strengths is desirable JJi]. We nevertheless believe that 
our results are useful because strong vortex tubes play an 
important role in small-scale intermittency. Their role in 
energy dissipation is also expected to be important. 
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